Nature of bonding and electronic structure in MgB2, a boron intercalation 

superconductor 
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Chemical bonding and electronic structure of MgB2 , a boron-based newly discovered superconduc- 
tor, is studied using self-consistent band structure techniques. Analysis of the transformation of the 
band structure for the hypothetical series of graphite - primitive graphite - primitive graphite-like 
boron - intercalated boron, shows that the band structure of MgB2 is graphite-like, with n bands 
falling deeper than in ordinary graphite. These bands possess a typically delocalized and metallic, as 
opposed to covalent, character. The in-plane a bands retain their 2D covalent character, but exhibit 
a metallic hole-type conductivity. The coexistence of 2D covalent in-plane and 3D metallic-type 
interlayer conducting bands is a peculiar feature of MgB2. We analyze the 2D and 3D features of 
the band structure of MgB2 and related compounds, and their contributions to conductivity. 



The^ recent discovery of superconductivity (SC) in 
MgB2cLphas stimulated significant interest in this 
systemoLl. One of the first questions is whether MgB2 
represents a new class of superconductors, or it may be 
related to other known classes in terms of its bonding and 
electronic properties. The crystal structure of MgB2 ma*; 
be regarded as that of completely intercalated graphitecl 
with carbon replaced by boron (its neighbor in the peri- 
odic table). In addition, MgB2 is formally isoelectronic 
to graphite. Therefore, understanding the bonding and 
electronic properties of MgB2 and its relation to graphite 
intercalation compounds (GIG), some of which also ex- 
hibit SG, is desirable and is the subject of this paper. 

The search for HTSG in carbon phases started in 1965 
when Tc = 0.55 K in KGg was reported, and SG was 
subsequently explained in terms of the interaction of tt- 
derived bands rwith the longitudinal phonon modes of the 
graphite layercl. Similar conclusions for the electronic 
structure at the Fermi level were derived for another jGIG, 
LiGe. The highest T^ achieved for a GIG was 5 KB. A 
parallel development was that of SG in Bechgard salts, 
the organic charge-transfer systems. Within this family 
of organic materials exhibiting carbon 7r-band conductiv- 
ity, Tc has been raised to 12.5 kQ. It should be noted that 
the majority of these compounds are unstable at normal 
conditions and high-pressure experiments are most com- 
mon. This part of GIG studies was reviewed in Ref. ||. 

Later research shifted to the area of carbides and ni- 
trides of transition metals where metallicity and conduc- 
tivity are mostly due to transition metal atoms. Only in 
1991 was it shown thai alkali-doped Ggo exhibits SG with 
maximum Tc of 33 KB. It was the first representative of 
a 3D carbon-based metal with modified (compared to 
graphite) coupling of lattice vibrations to the conduction 
electrons due to the curvature of the molecule but still 
with TT-band derived states at the Fermi level. 

We have calculated the electronic structure of MgB2 
and some related compounds using the Stuttgart LMTO- 
TB (ASA) code and full-potential LMTO (FLMTO) 



methodcl. A local exchange-correlation potential has 
been used in both cases (addition of gradient correc- 
tions did not change the results). It appears that a 
general analysis of energy bands may be safely done in 
ASA, while reliable treatment of charge densities and 
anisotropy of transport properties of these compounds 
requires more accurate FLMTO calculations. 

In order to examine the relation between the band 
structure of MgB2 and that of graphite we studied the 
following hypothetical sequence of intermediate materi- 
als: carbon in the 'primitive graphite' (PG) lattice with 
no displacement between layers as in MgB2 and in GIGs, 
using graphite lattice parameters; boron in the PG lat- 
tice with a as in MgB2 and c/a as in graphite; boron 
in the PG lattice with a and c/a as in MgB2; LiB2 in 
the same structure; MgB2 itself. The results of some of 
these calculations are shown in Fig. 1. The band struc- 
ture of PG carbon shown in Fig. la is very similar to that 
of graphite (see e.g. Ref. |l^ and references therein) with 
the appropriate zone-folding for a smaller unit cell. (This 
is quite natural because of the weak interlayer interac- 
tion.) Boron in the same lattice but scaled to match the 
MgB2 in-plane lattice parameter (not shown) has nearly 
identical bands with the energies scaled by the inverse 
square of the lattice parameter, in^agreement with com- 
mon tight-binding considerationald. Fig. lb shows the 
natural enhancement of the out-of-plane dispersion of the 
TT bands when the interlayer distance is reduced. Figs. 
Ic and Id demonstrate that 'intercalation' of boron by 
Li or Mg produces a significant distortion of the band 
structure, so that the role of the intercalant is not simply 
one of donating electrons to boron's bands (which would 
return the band structure to that of PG carbon shown 
in Fig. la). The main change upon intercalation is the 
downward shift of the tt bands compared to a bands. For 
Li this shift of ~ 1.5 eV is almost uniform throughout the 
Brillouin zone. Replacement of Li by Mg shifts the tt 
bands further, but this shift is strongly asymmetric in- 
creasing from ~0.6 eV at the F point to ~2.6 eV at the 
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A point. In addition, the out-of-plane dispersion of the 
tJ bands is also significantly enhanced. In LiB2 the filling 
of the antibonding a band is nearly the same as in PG 
boron, while in MgB2 the Fermi level shifts closer to the 
top of the antibonding a band. 

The Fermi sm^face (FS) for MgB2 obtained in the 
FLMTO method is practically identical to Fig. 3 of 
Ref. ^. The incompletely filled antibonding a band with 
predominantly boron p^y character forms two hole-type 
cylinders around the F — A line, while the tt bands form 
two planar honeycomb tubular networks: one electron- 
type one centered at fc^ = and another similar, but 
more compact, hole-type one centered at — ir/c. The 
lowering of the tt band in MgB2 compared to PG boron 
is connected with charge accumulation in the interstitial 
region. This lowering is greater at the A point compared 
to the F point, because the interstitial charges favor an- 
tisymmetric overlap of the boron's pz tails. 

The nature of bonding in MgB2 may be understood 
from the charge density (CD) plots shown in Fig. ^j. As 
seen in Fig. Ha, bonding in the boron layer is typically 
covalent. The CD of the boron atom is strongly aspher- 
ical, and the directional bonds with high CD are clearly 



seen (see also Ref. 
distribution in the 



121). It is worth noting that the CD 
Doron layer is very similar to that in 
the carbon layer of graphite (see e.g. Ref. [l^). This di- 
rectional in-plane bonding is also obvious from Fig. ^ 
showing the CD in a cross section containing both Mg 
and B atoms. However, Fig. |2|d also shows that a large 
amount of valence charge does not participate in any co- 
valent bonding, but rather is distributed more or less 
homogeneously over the whole crystal. Further, Fig. ^jc 
shows the difference of the CD of MgB2 and that of hy- 
pothetical NaB2 in exactly the same lattice. Not only 
does it show that one extra valence electron is not ab- 
sorbed by boron atoms but is rather delocalized in the 
interstitials; it also shows that some charge moves out- 
ward from boron atoms and covalent in-plane B-B bonds. 
Fig. ||d shows the CD difference between the isoelectronic 
compounds MgB2 and PG carbon (C2). In MgB2, the 
electrons see approximately the same external potential 
as in C2, except that one proton is pulled from each C 
nucleus and put at the Mg site. It is evident that the 
change C2— >MgB2 weakens the two-center a bonds (the 
charge between the atoms is depleted) and redistributes 
it into a delocalized, metallic density. 

Thus, the structure of MgB2 is held together by 
strongly covalent bonding in boron sheets and by delo- 
calized, 'metallic-type' bonding between these sheets. A 
peculiar feature of this compound is that electrons partic- 
ipating in both of these bond types provide comparative 
contributions to the density of states (DOS) at the Fermi 
level (see below). This distinguishes MgB2 from closely 
related GICs where covalent bonds in the carbon sheets 
are always completely filled, while the nearly cylindrical 
parts of the FS commonly found in those compounds are 
formed by carbon-derived tt bands which su-e also respon- 
sible for conductivity in pristine graphiteU. 



Because of the coexistence of two different types of 
chemical bonds, it is desirable to find the contributions 
to the total DOS and transport properties from sepa- 
rate pieces of the FS originating from 2D covalent and 
3D metallic-type bonding. Such decomposition is shown 
in Fig. H for the total DOS and for the in-plane {xx) 
and out-of-plane {zz) components of the tensor Uafj = 
^ VaVp5{£ — Ep)<fi\i., with Va being a component of 
the electronic velocity. One can see that the 3D (metallic- 
type bonding) and cylindrical (covalent bonding) parts 
of the FS contribute, respectively, about 55% and 45% 
to iV(0). N{e) for the hole-type zones rapidly decreases 
with increasing e and already at £ ~ 0.8 eV above Ep 
the total DOS is almost completely determined by the 
3D electron-type band. The latter contribution is almost 
constant and probably is not related to the change of 
SC properties under pressure or with doping. The cor- 
responding contribution to a exceeds all other contribu- 
tions (more than 50% for azz and axx) and is virtually 
isotropic. Holes in the in-plane B-B covalent bands (two 
cylinders at the FS) as expected have clearly anisotropic 
behavior contributing about 30% to axx and virtually 
nothing to CTzz- The 3D hole- type part of the FS is also 
notably anisotropic with predominantly z-axis conduc- 
tivity. The total a has a relatively small anisotropy at 
Ep with axx = 2.66 and azz — 2.18 Ry • (a/27r)^. 

The main challenge is to suggest another possible SC 
with a higher Tc. We studied several similar systems 
which, as we believe, may be interesting for experimen- 
tal studies. Several isoelectronic compounds demonstrate 
very different behavior. Hypothetical BeB2 and ZnB2 in 
MgB2 structure are very similar to MgB2 in terms of 
the electronic properties. Our ASA calculations of these 
compounds with lattice parameters of MgB2 produced 
very similar band structures with nearly identical FS. 
Presumably a smaller radius of the Be ion and a larger 
one of the Zn ion will generate some difference but we 
expect that these compounds (if existing) would have 
properties similar to MgB2. Anpther isoelectronic sys- 
tem is lithium borocarbide LiBCtJ, which, according to 
our ASA results, is a perfect insulator, so that any sub- 
stitution of C by B will lead to metallic behavior. (In 
general, the obtained band structure of this system is 
very similar to that of BNIiJ.) Experimentally very small 
conductivity was observed in LiBCllj. 

Addition of one electron to MgB2 corresponds to AIB2. 
This compound also has metallic conductivity but its FS 
does not have any cylindrical partsi-as in MgB2. Exper- 
imentally no Tc has been observcdila which makes elec- 
tron doping an unlikely mechanism to increase Tc. Un- 
der such circumstances it would be most logical to try 
hole doping by replacing Mg by Na or Li. Unfortunately 
such compounds are unstable. We ace familiar only with 
one experimental report that NaC2Lj is a weak SC under 
pressure and is very unstable. However, it may be possi- 
ble to form an alloy (Mg,Na)B2 with modest amounts of 
Na. It is evident from Fig. |^ that the hole-type parts of 
the FS change dramatically with electron filling. 
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In summary, MgB2 represents the first relatively high 
Tc superconducting compound of boron and has very 
peculiar bonding characteristics compared to GIC. We 
showed that the electronic structure of MgB2 has both 
similarity with and notable differences to the GICs. 
These features provide a basis both for further studies 
of normal and SC states of MgB2 and for engineering 
new SC compounds. 
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FIG. 1. Band structures of: (a) top left: primitive (AA 
stacking) graphite (PG), a = 2.456A, c/a = 1.363; (b) top 
right: PG boron, a = 3.085A, c/a = 1.142 (as in MgB2); 

(c) bottom left: LiB2 in MgB2 structure, same a and c/a; 

(d) bottom right: MgB2, same a and c/a. Energy is in eV 
relative to Ef- The order of occupied bands in the P point is: 
a bonding, tt bonding, a antibonding (double degenerate). 



FIG. 2. Pseudocharge density contours obtained in 
FLMTO. The unit cell is everywhere that of MgB2. Dark- 
ness of lines increases with density, (a) MgB2 in (0002) plane 
passing through B nuclei; (b) MgB2 in (1000) plane pass- 
ing through Mg nuclei at each corner of the figure. B nu- 
clei occupy positions (1/3,1/2) and (2/3,1/2) in the plane of 
the figure. The integrated charge of the unit cell is 8. (c) 
(1000) plane, difference in smoothed density, MgB2 minus 
NaB2. The integrated charge of the unit cell is 1. (d) (1000) 
plane, difference in smoothed density, MgB2 minus PG car- 
bon. The integrated charge of the unit cell is 0. In (c) and 
(d), dotted lines show negative values. 
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FIG. 3. Rigid band results for MgB2: (a) Total DOS in 
eV~^/celI; (b) a^x and (c) CTzz in Ry • (a/27r)^ (for definition 
of act0 see text), with contributions from different parts of the 
FS. Thick solid lines: total; thin dashed: 3D electronic part; 
dotted: 3D hole part; thick dashed: internal cylinder; thin 
solid: external cylinder. Energy is in eV relative to Ef- 
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